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ABSTRACT 

The large number of observed exoplanets (> 700) provides important constraints on their origin as 
deduced from the mass-period diagram of planets. The most surprising features in the diagram are 
1) the (apparent) pile up of gas giants at a period of ^ 500 days 1 AU) and 2) the so-called 
mass-period relation which indicates that planetary mass is an increasing function of orbital period. 
We construct the evolutionary tracks of growing planets at planet traps in evolving protoplanetary 
disks and show that they provide a good physical understanding of how these observational properties 
arise. The fundamental feature of our model is that inhomogeneities in protoplanetary disks give rise to 
multiple (up to 3) trapping sites for rapid (type I) planetary migration of planetary cores. The viscous 
evolution of disks results in the slow radial movement of the traps and their cores from large to small 
orbital periods. In our model, the slow inward motion of planet traps is coupled with the standard 
core accretion scenario for planetary growth. As planets grow, type II migration takes over. Planet 
growth and radial movement are ultimately stalled by the dispersal of gas disks via photoevaporation. 
Our model makes a number of important predictions: that distinct sub-populations of planets that 
reflect the properties of planet traps where they have grown result in the mass-period relation; that 
the presence of these sub-populations naturally explains a pile-up of planets at 1 AU; and that 
evolutionary tracks from the ice line do put planets at short periods and fill an earlier claimed "planet 
desert" - sparse population of planets in the mass-semi-major axis diagram. 

Keywords: accretion, accretion disks — turbulence — Methods: analytical — planets and satellites: 
formation — protoplanetary disks — Planet-disk interactions 



1. INTRODUCTION 

The availability of large samples of exoplanets is being 
used to const rain theories of plane t formation in a sta- 
tistical sense (jUdrv fc Santosl[2007[ ). The standard theo- 
retical tools for this are the so-called population synthe- 
sis mod els (Ida & Lin 2004, 2008a; M ordasini et al.il2009l: 
llda fc Linll201Q ). wherein gas giants are considered to be 
formed by two main suc cessive processes: the formation 
of cores by runaway fe .g.lWetherill fc Stew art 1989) and 
oligarchic growth (e.g. iKokubo fc IdalliggS). followed b y 
gas accretion onto the cores fe.g. IPollack eFallliggel) . 
This mode of forming gas giants is referred to as the 
core-accretion scenario. The orbits of these accreting 
protoplanets are regulated by planetary migration that 
eventually d etermines the radial distribution of planets 
(|WardllT997f ). The spirit of population synthesis mod- 
els is to hypothesize that the diversity in the properties 
of observed exoplanets reflects the range of the (initial) 
disk environments in which planets are born. Fine tun- 
ing of the efficiency of various physical processes such as 
migration rates allows one to qualitatively reproduce the 
observations summarized in the mass-period diagram. 

Despite the success of these models, a single complete 
theory of planet formation that can reproduce the archi- 
tecture of any (exo)planetary system including our Solar 
system is still unknown. In particular, it is unclear as to 
the physical origin of several key observations: the (ap- 
parent) pile-up of planets at 1 AU and the mass-period 
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relation which shows that planetary mass increases with 
period (see Fig. [1J0 Furthermore, there is a signifi- 
cant discrepancy between the theories and observations: 
the recent population synthesis models claimed that a 
planet desert - a region in the mass-period diagram with 
a lower population of exoplanets - is present in the range 
of planetary mass {5Mq < Mp < 50 Mm) and of their 
semi-major axis (0.04 AU< r < 0.5 AU) (llda fc Linll2004 
'2008b'13 whereas many exoplanets are already observed 
there (see the black rectangle in Fig. [T|). 

In this paper, we address how inhomogeneities in pro- 
toplanetary disks can account for the observed trends. 
For this purpose, we constructed and followed evolution- 
ary tracks of planets that grow at disk inhomogeneities. 
More specifically, we compute planetary growth and mi- 
gration in protoplanetary disks that evolve with time due 
to disk viscosity and photoevaporation of gas, by track- 
ing the movement of disk inhomogeneities such as dead 
zones, ice lines and heat transitions. The fundamental 
contribution of disk inhomogeneities to theories of planet 

^ Observations prefer orbital periods while semi-major axes are 
more natural in theoretical calculations. Since they are translatable 
through some analytical relations, we converted the observational 
data of Mavor et al. (2011) from periods to semi-major axes using 
their published data of periods, planetary mass, eccentricities, and 
the amplitude of the radial velocities. Thus, we mainly use semi- 
major axes rather than periods. 

^ Recently, Ida & Lin| 112010 ') succeeded in reproducing the pop- 
ulation of low mass planets with short orbital radii by adding an- 
other physical process - mergers of protoplanets - that takes place 
'caj.fter the gas disks are severely depleted. As shown below, on the 
contrary, our model is able to explain the population within the 
same framework of forming gas giants. 
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Figure 1. Observed exoplanets using the radial velocity tech- 
nique (denoted by black pluses). The data are obtained by the 
CORALIE and HARPS surveys, both of which are carried out 
through modest- and high-resolution spectrographs at La the Silla 
Observatory in Chile. We took the data from Mayor et al. ( 2011), 
wherein ~ 150 observed exoplanets are selected from larger samples 
using consistent observational criteria for their statistical analyses. 
Thus, the data are well defined in order to discuss the statistical 
properties of exoplanets. Also, MpSini is plotted, since the incli- 
nations i are unknown. We have converted the data from orbital 
periods to semi-major axes for direct comparisons with our results 
(see Fig. \5^. The host stars are F, G, or K stars. Jupiter (the red 
cross) and Saturn (the red star) are also shown for the reference. 
The thick black lines denote the amplitude of redial velocities of 
1 m s~^ and 10 cm s ~^ . The amplitude of 10 cm s ~^ is not 
achieved yet even in the HARPS survey while 1 m s~^ is well in 
hand. The data show the two trends; the (apparent) pile up of 
gas giants at ~ 1 AU and the mass-period relation wherein plane- 
tary mass is an increasing function of period (particularly beyond 
1 AU). Earlier papers (Ida k Lin 2004, 2008b) predicted a planet 
desert demarcated by the black rectangle that covers in the range of 
planetary mass (SM^ < Mp < 50Mg)) and of their semi-major axis 
(0.04 AU< r < 0.5 AU) in the diagram. The recent observations 
populate the desert. 



formation here is that they give rise to trapping sites 
for rapid type I planetary migration of cores of gas gi- 
ants (often referred to as planet trap s in the hterature, 
Masset et all 120061: llda fc LinI l2008bl: iMatsumura et alj 
2009; Has egawa fc Pudritzll2010aHLvra et al.ll2010[) . Fol- 
lowing the viscous evolution of disks, planet traps gradu- 
ally move inwards in their disks, taking the trapped cores 
with them. The trapping and transport of cores is a cen- 
tral feature of our models in which protoplanets accrete 
gas as they move with the traps. We will show below 
that a semi-analytical model, wherein these two effects 
of planet traps, further planetary growth and subsequent 
type II migration that is terminated by photoevapora- 
tion of gas disks are all combined, can provide natural 
explanations of a number of the important observational 
properties: 1) the origin of the observed mass-period re- 
lation, 2) the origin of the pile up of observed gas giants 
at ^ 1 AU, 3) the origin of low-mass planets distributing 
in the earlier claimed planet desert, 4) prediction of a 
new planet desert that originates from different physical 
processes than the earlier desert. 

The plan of this paper is the following. In §[51 we sum- 
marize under what conditions planet traps are generated 
and which tidal torque and disk property play the most 



crucial role for activating the traps. In § [3J we describe 
disk models that are used for specifying the properties of 
disk inhomogeneities while, in § |4l we discuss how these 
inhomogeneities evolve with time following viscous evo- 
lution of disks with photoevaporation of gas. In § [5l we 
derive the characteristic masses of planets that are cap- 
tured at planet traps and how these masses define the 
mode of planetary migration. In § [6l we synthesize the 
above treatments and develop a semi-analytical model of 
planetary growth and migration affected by planet traps 
for constructing evolutionary tracks of growing planets 
in the mass-period diagram. We present our results and 
compare them with the observations in § [T] The general 
reader may wish to skip to § [7] for a non-technical, astro- 
physical discussion of the results. Parameter studies are 
performed in § [51 § [9] is devoted to our discussion and 
conclusions. 

2. ORIGINS OF PLANET TRAPS 



Planet traps, a term first coined by iMasset et alJ 
(2006), are one of the keys to resolving the long-standing 
problem of rapid type I migration that can lead to the 
loss of any planetary system to the host stars within 
~ 10^ years. This arises due to the high efficiency 
of angular momentum transfer between (proto)p lanets 
and their natal disks (|Wardlll997t iTanaka et aLr r2002). 
The basic idea of planet traps lies in the fact that 
the direction of migration can switch from inwards 
to outwards when planets migrate through the disks 
that have some kind of inhomogeneities. This is 
the combined consequence of the high sensitivity of 
type I migration to disk properties (iTanaka et alJ 120021 : 
iPaardekooper et al.ll2010HHasegawa fc Pudritdl2011al lbl) 
and the density and temperature modification s pro- 
duced by the disk inhomogeneities (iD'AlessicTet al.l 
nigg: iMenou fc Goo dman"2004l: IMatsumura et al.ll2007l: 
IHasegawa fc Pudritz 2010b) . Since a complete discus- 
sion of how disk inhomogen eities give rise to planet traps 
is presented elsewhere (e.g. IHasegawa" fc Pudritzl[20TTbl 
hereafter Paper I), we simply summarize the typical disk 
configurations with which planet traps are created and 
the responsible tidal torques and disk properties by which 
planet traps are activated with such disk configurations 
in Table [H In this paper, we focus on dead zones, ice 
lines, and heat transitions which all become planet traps 
(see §[3|for their definitions, also see Paper I). 

3. DISK MODELS 

A significant number of analytical and numerical 
studies have shown that "realistic" disks a re likely 
to possess several kinds of inhomogeneities (Gammid 
1996; D'AlcssioeLaLl 119981 : iMenou fc Goodm an 200i 
(Min et al.ll201li r which can activate planet traps. We 
first briefly describe our disk models that serve as the 
basis for specifying the properties of the disk inhomo- 
geneities such as their positions and surface densities. 
We refer the reader to Paper I for the complete discus- 
sion. 

We adopt the standard models of steady accretion 
disks that have accretion rates modeled as 



(1) 



where Sn 



aCsH, Cs, and H are the surface 



density, the viscosity, the sound speed, and the pres- 
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Table 1 

Typical disk structures for generating planet traps 



Torque 



Disk structures 



Relevant disk properties 



Lindblad torque 

Vortensity-related corotation torque 
Entropy-related corotation torque 



t > 1.5 with s = 
s < —2 with t = 

s > 1 with t ^ 
t < —IT with s 



— 1 Dead zones 

—0.5 Dead zones and ice lines 
-0.5 Inner edge of disks 

= — 1 Viscous heating 



Simple power-law disk structures are assumed, that is, the surface density of gas S oc r° and the disk temperature 
T oc r*. We refer the reader to Paper I for a more complete discussion. 



sure scale height of gas disks, respectively. The fa- 
mous a— prescription is ass umed for characterizing th e 
strength of disk turbulence (jShakura fc Sunvaevlll973l) . 

3.1. Positions of disk inhomog entities 

Adopting the standard disk model, we can estimate the 
positions of disk inhomogeneities. These positions are 
crucial because (proto)planets that undergo rapid type I 
migration will get trapped there. We simply summarize 
the positions here and refer the reader to Paper I for the 
complete derivations (also see Table [J). 

There are generally of three kinds of disk inhomo- 
geneities: dead zones, ice lines and heat transitions 
(Paper I). Dead zones are present in the inner re- 
gion of disks where high energy photons such as X- 
rays from the central sters and cosmic rays cannot 
penetrate (iGanimiel 119961 : iMatsumura fc Pudritd 120061 : 
lllgner fc NelsonI 120061 )" The defining feature of the dead 
zones is the low amplitude of turbulence there that re- 
sults from the poor coupling of the magnetic field with 
weakly ionized disks (so that magnetorotational insta- 
bilities (MRIs) are suppressed there, e.g. Balbus 200»|). 
Ice lines at which the disk temperature is low enough to 
trigger condensation of molecules such as water are the 
most famous and indispensable of disk inhomog eneities 
(|,Tang-Conden fe Sasselovl[2nnl IMin et al.ll20Tlh . They 
play an irnporta n t role in population s ynthesis models 
(|Ida fc LinI 120041 : iMordasini et al.ll2009t ). Heat transi- 
tions are also well recognized in the literature and arise at 
that radius at which the main heat sour ce changes from 
viscous heating to stellar irradiation ([D'Alessio et al.l 
119981 iMenou fc GoodmanI [2004. Paper I). We have re- 
cently demonstrated that the heat transition becomes a 
pla net trap based on an alytical arguments (Paper I, also 
see lKretke fc Linll2012|) . The validity of the heat transi- 
tion traps h as been recently confirm ed by hydrodynamic 
simulations ()Yamada fc Inaball2012[ ). 

In principle, the structure of dead zones can be spcci- 
fied by solving the ioniz ation eq uations ( Sano_ et al 200Q; 
IMatsumura fc Pudritd 120061 : lllgner fc' Nelso3 (200l l. 
Nonetheless, the resultant structures depend sensitively 
on disk parameters that are difficult to determine 
through the observations. Therefore, we adopt a param- 
eterized treatment of dead zones (Kretke fc Lin 2007; 
llda fc Lid l2008al : IMatsumura et al.n 2009. Paper I) in 
which the effective a in the layered region can be given 
as 



^AaA + (Sg - TjA)aD 



(2) 



where a a and ao are the strength of turbulence in the 
active and dead layers, respectively, and the surface den- 



sity of the active layer is modeled as 

SA 

I I- \ 

— ^Aofice 



(3) 



where = 1 AU is the characteristic radius, fice can 
contain the effects of ice lines. More specifically, fice 
represents possible reductions of at the ice line that 
originate from the complex interplay between ice-coated, 
sticky du st grains and the absor ption of free electrons 
by them (jSaiio et al.ll2000Hlda fc Lin 2008li Paper I, see 
a more complete discussion). Thus, the structures of 
dead zones are controlled only by 'Sao and sa in this 
formalism. This approach is useful because it enables one 
to investigate how important the structure of dead zones 
is for understanding the population of planets by simply 
varying parameters, T,ao and sa (see § [5]). Assuming 
stationary disk models (see equation (|T])), the surface 
density of gas is given as 



o^A - ap 
CUD 



(4) 



With equation Q in hand, the position of a dead zone 
trap is given as 



rdz 
To 



M 



(5) 



where Hq and ilo are the pressure scale height and Kep- 
lerian frequency at r — tq. This can be derived from the 
assumption that the outer edge of dead zones is specified 
around 'Sa Sg/2- 
The position of the ice line of molecular species k is 




2/9 



(6) 

where Tm^k is the disk midplane temperature below 
which molecules k can condense, So = 2 x 10^^ is the 
opacity at the ice line of the molecule, /ig is the mean 
molecular weight of the gas, fc^ is the Boltzmann con- 
stant, and 7 = 1.4 is the adiabatic index. This is given 
by the recent results which show that the viscous heat- 
ing (rather than stellar irradiation) is gener ally domi- 
nant for determining the position of ice lines (jMin et al.l 
I2011L Paper I). This expression is applicable for any 
molecules. Nonetheless, we focus on water ice lines here, 
since they are likely to be the most important molecule 
for understanding the observed mass-period relation (Pa- 
per I). For ice lines of water, T^-^ o(^^') ~ ^'^^ ^ 
(|Jang-Condell fc Sasselovll200l . 
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Figure 2. Time evolution of the positions of three disk inhomo- 
geneities for a disk around a classical T Tauri star (CTTS). The 
dead zone is denoted by red, the ice line by green, and the heat tran- 
sition by light blue. Following the viscous evolution of disks, these 
zones move inwards at different rates. This results in the complex 
behaviors; convergence, disappearance, and re-appearance. Their 
end points are determined when photoevaporation of the gas disk 
takes place. 



For the case that ice hnes are located within dead 
zones, the position of the trap needs to satisfy the fol- 
lowing condition: 



> h{rdz) 

Tdz \ aA- CUD , 



(7) 



Finally, the position of heat transitions is written as 



1 / ro 



3/7 




14/15 

(8) 



OCM28/45, 



where kq = 2 x 10 ^ is the opacity at the heat transitions, 
i?* is stellar radius. 



T, 



mO 



2/7 



GAL 



1/7 



knR* 



(9) 



(10) 



and are stellar effective temperature and mass, 
respectively, and G is the gravita tional constant. We 
have a dopted analytical models of iChiang fc GoldreichI 
(|1997| ) for the temperature of the disk midplane heated 
by stellar irradiation. Planet traps arising from the heat 
transitions are active only if r^t > rdz- 

By comparing the positions of each disk inhomogene- 
ity (see equations ([5]), ([6]), and one immediately 

observes that disk evolution, which lowers the accretion 
rate M, moves them inwards, but at different rates (see 
Fig. [2]). This is important for understanding the ob- 
served mass-period relation. 



3.2. Characteristic surface densities at disk 
inhomogeneities 

We estimate the characteristic surface densities at 
the disk inhomogeneities, following the above formula- 
tion. Although the detailed structures of disk inho- 
mogeneities remain to be simulated, a number of an- 
alytical and numerical studies based on the standard 
viscous di sk theory cla rified their chara cteristic struc- 
ture (Menou fc Goodman 2004; Matsu mura et all 120071 : 
llda fc Linll2008bD . These surface densities are utilized 
for deriving the characteristic masses of planets and fol- 
lowing evolutionary tracks of planets that grow in planet 
traps. 

As briefiy mentioned above, the outer edge of dead 
zones is determined by T,a ^3/2- Substituting this 
condition into equation (|4]), we find that the characteris- 
tic surface density at rdz is given as 



2M 



-•g.dz 



37r{aA + aD)r'^h'^^ 



(11) 



where h = H/r is the aspect ratio. 
At the ice lines, the surface density is approximately 

written as 

~ STTaor^h^n- ^^^^ 

We took the mean value of a as ~ au. As mentioned 
before, this assumpt ion is based on the recent extensive 
studies of ice lines (ISano et al.l[2000l: llda fc Linll2008bL 
Paper I). These studies indicate that ice lines can be 
regarded as a localized dead zone. 

On the other hand, the magnitude of turbulence at the 
heat transition is expected to be high enough to assume 
that disks are fully turbulent. As a result, the charac- 
teristic surface density at the heat transitions is given 
as 

Sna^r^h^n- ^^^^ 

4. TIME EVOLUTION OF DISKS AND THEIR 
INHOMOGENEITIES 

Time evolution of protoplanetary disks is established 
by the combination of viscous turbulence and the photo- 
evaporation of gas. These agents regulate the movement 
of disk inhomogeneities. We present our treatments of 
them. 

4.1. Viscous evolution 

Viscous turbulence is the dominant driver of disk evo- 
lution (e.g. Armitagc 2011). We adopt similarity solu- 
tions (|Lvnden-Bell fc Pringldll974( ) for constraining the 
relation between the accretion rate M and time r. Sim- 
ilarity solutions are derived from the conservation of 
angular momentum of disks (also see iHartmann et al] 
[199 81. Considering a disk that has mass Md and a char- 
acteristic disk radius Rc, its angular momentum Jd can 
be written as 

Jd « MdR'J^ (14) 

Following time evolution where disk material is accreted 
onto the central star, equation p4)) ensures that Rc is an 
increasing function of time (since Jd is roughly constant 
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and Md steadily decreases with time). In a simplified 
analysis, the expansion rate of Rc can be written as 

(15) 



dt 



Rc 



where t^is — /{iv) is the viscous timescale. Assuming 
a power-law structure for the disk temperature (T oc r*), 
equation ([T5l) gives 



Rd oc T<l/2-t) , 



(16) 



and the total disk mass Md decreases as (see equation 

m) 

As a result, the accretion rate is related to time through 
the following relation; 



A I oc r '^TT^ . 



(18) 

Combining the observations which show that the me- 
dian accretion rate for classical T Tauri stars (CTTSs) of 
age ~ 1 Myrs is ~ IQ~^Mq yr'^ fHartm ann et al.lll998l) 

and that M oc (jCalvet et al., .2004: .Muzerolle et alj 
[2005I) . we have the following scaling law for accretion 
rates; 



106 yr 



0.5M, 



(19) 

where we have assumed that the typical mass of CTTSs 
is ^ 0.5Mq and introduced a dimensionless factor face- 
This factor can be utilized for varying M and investigat- 
ing the subsequent consequences on disk evolution and 
planet formation. 

4.2. Photoevaporation 

It is still unclear how gas dis ks disperse in the fi- 
nal stages of their evolution (e.g. IArmitagil2011[ ). One 
of the leading mechanisms is photoevaporation which 
arises from heating up gas by high energy photons from 
the surrounding stars and sub sequent evaporation of 
gas due to the therma l pressure (jHollenbach et al.lll994l : 
iJohnstone et al.l [19981 ). In principle, photoevaporation 
rates are determined by the complex interplay between 
physical and chemical processes that take place in pro- 
toplanetary disks being irradiated by their central and 
nearby massive stars from fa r-UV (FUV) to extreme- 
UV (EUV) and up to X-rays (jCorti HollenbachllMil 
references herein) . Recent extensive studies have investi- 
gated how effective photoevaporation is in the dispersal 
of gas disks. In our models, we adopt a simple scaling 
law to represent the effects of photoevap orat i on. 

Following the treatment of Adams et al.l ()2004 see 
their Appendix for the complete derivation), photoevap- 
oration rates can be scaled as 

Mpe = fpeNctigC.Tg (^) CXp (-g) , (20) 

where fpe is a dimensionless factor of order unity, Nc is 
the critical column density of gas that is heated by stellar 
radiation, and the gravitational radius Vg is given as 



GM, 



knT 



100 AU 



/ T 



V 1000 K 



Ah 

lAIr. 



(21) 



Utilizing some of the most advanced results of pho- 
toevap oration, we further simplif ied equation pO)). Re- 
cently, iGorti fc HoUenbachl (|2009f) have investigated pho- 
toevaporation of gas disks by taking into account radia- 
tion of a central star that covers FUV, EUV and X-rays, 
and found that FUV heating plays the dominant role for 
inducing photoevaporation at r > 3 AU. This can be 
understood by the fact that photoevaporation rates are 
determined by the product of the gas temperature and 
density. EUV heating that leads to ionizing atomic hy- 
drogen results in higher gas temperatures (^ 10* K) than 
FUV heating (-- 10^-10^ K). Nonetheless, the ioniza- 
tion front above which EUV heating dominates can only 
penetrate the disk atmosphere where gas density is much 
lower than that where FUV heating becomes dominant. 
As a result, FUV-induced photoevaporation rates exceed 
EUV-induced ones. 

When photoevaporation is established mainly by FUV 
heating, the heated outgoing flow acts as an additional 
sourc e of opacity fo r the F UV photons (jJohnstone et al.l 
Il998t lAdams et al.l [200l . Consequently, the critical 
column density heated up by FUV satisfies the self- 
regulation relation; 



TPUV = (^FUvNc ^ 1, 



(22) 



where (Jfuv ~ 8 x 10^^^ cm^ is the reasonable cross sec- 
tion of dust grains for the FUV photons. This enables us 
to specify Nc in equation ([20]) . Also, the peak of FUV- 
induced photoevaporation rates is attained around 0.1- 
0.2rg, rather th an Tq that is valid for EUV-induced pho- 
toeva poration (|Adams et al.l 120041 : IGorti k, HoUenbachl 
l200a . This can be again explained by the combination 
of the gas temperature and density, and is also confirmed 
by equation (^0)1 . 

Collecting the above arguments, we obtain a simplified, 
but physically motivated scaling law for photoevapora- 
tion rates; 



Mpe ^ 2.3 X lO-^Mo yr- 



3 km s ^ 



Al, 
IMq, 
(23) 

and set that 



where we have used equations (|2ip and 
r = 0.1 rg in equation (|20|) . 

We note that equation ([23]) can be applied for pho- 
toevaporation rates induced by both EUV and X-rays 
despite the fact that it is derived from the physical con- 
sideration based on FUV radiation. This can be done 
by adjusting the dimensionless factor fpe that is deter- 
mined by the compariso n with more detailed simula tions. 
In fact, the conclusion of lGorti fc HoUenbachl (j2009f ) that 
FUV is the dominant source of photoevaporation is still 
a matter of debate in the literature. This is partly be- 
cause they relied exclusively on hydrostatic solutions for 
quantifying winds driven by FUV radiation (although 
hydiodynamical models are needed for precisely estimat- 
ing the winds), and partly because they adopted energy 
spectra whi ch even tually reduce the effects of X-rays. As 
shown by, Owen et al. (2010, 2011, 2012), photoevapora- 
tion rates induced by X-rays can attain ^ 10~^ M0 yr~^ 
for the most luminous X-ray sources. This high value 
is derived from employing observed Chandra spectra of 
TTSs and is comparable to the photoevaporation rate 
induced by FUV radiation. As a result, we intentionally 
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avoid specifying the dominant source of photoevapora- 
tion of gas. Instead, we consider the general effects of 
photoevaporation on planet formation by treating /pe as 
a free parameter. 

4.3. Photoevaporation of viscous disks 

We are now in the position to discuss the complete 
treatment of disk evolution. We assume that the accre- 
tion rate through the disk is constant in space and reg- 
ulated in time by equation ()19|) . As time goes on, disk 
material accretes onto the central star and the accretion 
rate decreases. This change in M drives the movement of 
the disk inhomogeneities (see equations ([5]), dS]), and dS])). 
When M becomes equal to the photoevaporation rate 
Afpe, represented by equation (|23|) . we assume the gas 
disks to disperse completely. Although this treatment is 
somewhat idealized, it can ac count for the more detailed 
simulations. As an example, iGorti et al.l ()2009l ) investi- 
gated the evolution of viscous protoplanetary disks that 
are photoevaporated by the FUV, EUV, and X-ray radi- 
ation from their central star. They showed that viscous 
turbulence controls the early stage of disk evolution and 
the total disk mass gradually decreases initially, which 
can be formulated by power-laws. Once the condition 
that M ~ Mpe is satisfied, the disk mass drops exponen- 
tially due to the combination of viscous evolution and 
photoevaporation of gas. They found that disks of initial 
mass O.IM© around ~ IM© have the lifetime of ^ 4 x 10^ 
years. 

We can derive similar results from our treatment by 
equating A'l with Mpf.'- the disk lifetime is approximately 
estimated as ~ 6 x 10^ years. Thus, our treatment is 
sufficient for the purpose of representing the evolution of 
protoplanetary disks that are regulated by both viscosity 
and photoevaporation. 

In summary, we reduce the surface density of gas and 
the accretion rates, following equation ([19)) . This results 
in the movement of the disk inhomogeneities. Also, we 
locate the final position of each disk inhomogeneity that 
is determined by the condition that M = Mpg. 

4.4. Parameters 

We summarize important parameters that establish 
the configuration and physical state of protoplanetary 
disks (see Table [2|). They are divided into three sets: 
stellar parameters (M*, T,, and i?*), the disk mass (Sao: 
sa, face, and t), and disk evolution (a^, cku, and fpe)- 
The disk evolution parameters can be translated into disk 
lifetimes. These three sets are fundamental to regulate 
planet formation in the ciisks, as confirme d in the popula- 
tion synthesis models (Idn \- Lin 2004; Mordasini et al.l 
|2009| ). We focus on disks around CTTSs and denote the 
set of the values given in Table [5] as our fiducial model. 
Adopting these parameters, disk evolution proceeds from 
Tint = 10^ year to the time at M = Mpe that defines the 
disk lifetime (rdisk)- A parameter study in which some 
of these quantities are changed is presented in § [S] 

4.5. Movement of planet traps 

We draw upon our comprehensive analytical study (Pa- 
per I) which showed that the gas surface density and 
temperature modifications induced by the disk inhomo- 
geneities are significant enough to reverse the direction 



of rapid type I migration. Therefore these positions are 
indeed trapping points of rapid type I migrators. 

Fig. [2] shows the movement of all the three disk in- 
homogeneities. As demonstrated by Paper I, they all 
move inwards, but at different rates. The inward move- 
ments arise from the viscous evolution and the resultant 
reduction of the surface density of the disks (see equation 
^9]) ). Also, different moving rates for the traps result in 
complex behaviors of the multiple inhomogeneities such 
as convergence (e.g. merging of the heat transition with 
the dead zone), disappearance, and re-emergence of in- 
homogeneities (e.g. the behavior of the ice line). When 
M equals Mpe, photoevaporation quickly disperses gas 
in the disk and hence the movement of the planet traps 
is terminated. This also determines the lifetime of the 
disk which is ^ 8.8 x 10^ years in this configuration. The 
behavior of the planet traps is crucial for understanding 
the observed mass-period relation later. 

5. CHARACTERISTIC MASSES 

We describe four characteristic masses that are impor- 
tant in our models (see Table [S]). Using the positions and 
characteristic surface densities at disk inhomogeneities 
given in § [3l these masses define the mode of planetary 
migration (see §[6|). Also, they result in a segment of the 
mass- semi- major axis diagram (see Appendix |A|) . 

5.1. Type I regime 

For planets of mass smaller than the gap-opening mass 
Mgap (see below), type I migration is the main agent that 
governs their orbital distribution. In our models, rapid 
type I migration is halted at the planet traps. Hence the 
location of type I migrators is predicted by the positions 
of disk inhomogeneities (Paper I, see Fig. [2]). 

It is important to define the minimum mass of planets 
that will be captured by the plan e t trap s. As demon- 
strated numerically by iLvra et all ()2010D , planets cap- 
tured in planet traps "drop-out" if the following condi- 
tion is satisfied: 

Tmig,I 



>1, 



(24) 



where Tmig.i is the timescale of type I planetary migra- 
tion and is the timescale that determines the moving 
rates of planet traps. This relation expresses the fact 
that trapped planets drop-out if the speed of type I mi- 
gration becomes less than that of the moving traps. In 
general, Tmig,i is scaled as 



MpTplp 

2T ' 



where 



fMp 



2 y 4q2 
lip 



(25) 



(26) 



with K-mig = 1 — 10, depending on the optical thickness 
of the disk (jPaardekooper et all I2O1O0 . For stationary 
accretion disk models (see equation (|I|)), the minimum 
mass of type I migrators that can be captured at planet 
traps is given as 



M, 



mig,! 



2K, 



(27) 
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Table 2 

Important disk quantities 



Symbols 


Meaning 


U i i OS 


M 


Stellar mass 


U.O IVIQ 




Stellar radius 


2.5 Rq 


T. 


Stellar effective temperature 


4000 K 




Surface density of active regions at r = ro 


20 g cm-2 




Power-law index of Sa(c< r"^ ) 


3 


t 


Power-law index of the disk temperature (T oc r*) 


-1/2 


face 


a dimensionless factor for M (see equation il9l) 


1 


OlA 


Strength of turbulence in the active zone 


10-3 


ao 


Strength of turbulence in the dead zone 


10-4 


fpe 


a dimensionless factor for Mpe (see equation II20II) 


1/3 



Table 3 

Characteristic masses 



Symbols 


Meaning 


Equation 


Mmig,I 


Minimum mass of type I migrators 


123 


^gap 


Gap-opening mass 


ll28t 


Merit 


Critical mass of type II migrators 


ll29t 


Mmax 


Maximum mass of planets 


|[30t 



We set = 10^ yrs, because the moving rates of planet 
traps are eventually regulated by disk lifetimes (see Fig. 
[2]) and the observations revealed that the disk lifetime of 
any CTTS disk is an order of Myrs. 

5.2. Gap-opening mass 

The gap-opening mass Mgap distinguishes type I mi- 
gration fro m type II and is well discussed in the U ter- 
ature (e.g. lWardlll997t iMatsumura k Pud ritz"2006'). It 
arises when a planet becomes sufficiently massive that 
the torque it exerts on the disk opens a gap. There are 
two main arguments for estimating Mgap. The first one is 
the Hill radius analysis: the Hill radius should be larger 
than the pressure scale height for maintaining gap forma- 
tion, otherwise gaps are closed by the gas pressure. The 
second argument arises from viscous disks. Disk viscos- 
ity that controls disk evolution plays the counteractive 
role for gap formation. Therefore, the tidal torque of a 
planet on their disks opens a gap if it exceeds the viscous 
torq ue. Summarizing these argum ents, Mgap is given as 
(e.g. IMatsumura fc Pudritzir2006[ ) 



gap 



3hl, Jmahl 



(28) 



5.3. Type II regime 

Planets of mass larger than Mgap open up a gap in 
their disks and undergo so-called type II migration. In 
the type II regime, we define two characteristic masses. 
One of them is the critical mass {Merit) above which the 
inertia of type II migrators is significant enough to pre- 
vent type II migration from proceeding as disks evolve 
(otherwise the timescale of type II migration is given as 
^ T,„:.s). This effect is also known as a damming 
effect (ISver fc ClarkellTQQaiWov et al.l [l999l The crit- 
ical mass Merit is defined by the local disk mass: 



Merit = ttE 



(29) 



where Sg.p is the surface density of gas disks at the po- 
sition of a planet (r = Vp). 



The other characteristic mass is the maximum mass 
of planets. In general, gas accretion onto cores of gas 
giants is no t fully terminated ev en if they form a gap in 
their disks ([Lissauer et al.|[2Q09l ). This suggests that the 
possible maximum mass of planets which start forming 
at time r can be estimated as 



Mrnax {t ) - 



drM 



= 5 X W^Hlefacc t 



(30) 



1 



M{Tdisk) 



W-^faccMQ yr- 



M(r) 



2 J \0.5Mq 

-1/(2(4-1)) 



-l/(2(t-l))- 



2t-l 

t-1 



lO-s/accMo yr- 



where equation ([T9| is used. 

In conclusion, the trapping regime is defined by Mgap 
and Mmigj in which type I migrators follow the move- 
ment of the planet traps while the type II regime is de- 
fined by Mmax and Mgap, wherein the radial distribution 
of planets is established by the type II migration (see the 
bottom panel of Fig. [TOl (Right) in Appendix [5| . 

6. EVOLUTIONARY TRACKS OF GROWING PLANETS IN 
PLANET TRAPS 

Armed with the positions and surface density of planet 
traps (§ [3]) and four characteristic masses (§ [5]), we 
now describe semi-analytical models of planetary growth 
and migration that are used for generating evolutionary 
tracks of accreting planets. 

6.1. Planetary growth 

The formation of gas giants is divided mainly into three 
stages ( Wctherill & Stewart 1989; Kokubo & Ida 1998|; 
IPollack et al.l 11996): formation of rocky cores through 
runaway and oligarchic growth (Stage I), the subsequent 
slow gas accretion of the cores and formation of envelopes 
surrounding them (Stage II), and collapse of the en- 
velopes and runaway gas accretion onto their cores (Stage 
HI). In order to model these three physical processes, we 
adopt the formulation of Ida & Lin (2004) who first at- 
tempted to understand the statistics of the observed ex- 
oplanets by carrying out population synthesis analyses. 
In this formulation, these processes are treated by sim- 
ple, analytical prescriptions that are derived from the de- 
tailed numerical simulations. In Appendix [B] we briefly 
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describe our treatments that slightly modify the origina l 
formulation, and refer the readers to llda fc LinI (|2004[ ) 
for a complete discussion. We also present a parameter 
study in Appendix [C] for confirming the validity of our 
tiny modifications. 

6.2. Orbital evolution of planets 

Orbital evolution of planets is governed by planetary 
migration that arises from tidal interactions ofthe plan- 
ets with the surrounding gaseous disks (IWMdl [19971: 
iTanaka et al.l |2002) . As discussed in § [5l the charac- 
teristic masses will classify planetary migration to four 
modes, depending on planetary mass: slower type I, 
trapped type I, the standard type II, and slower type 
II migration (see Table [3]). We discuss our treatments of 
them below. 

Slower type I migration: This mode is applicable if 
planetary mass is smaller than Mp < Mmigj (see equa- 
tion ([271) '). When planets satisfy this condition, the mi- 
gration rate of these planets is much smaller than the 
moving rate of gas that is regulated by disk viscosity. 
This is the reason why we call this mode of migration 
the slower type I migration. Therefore, we assume that 
these planets remain in the same position with time. 

Trapped Type I migration: When planets are in the 
trapping regimes, that is, Mmig.i < Mp < Mgap (see 
equations ([27| and (dH)), the radial positions of these 
planets follow the movement of planet traps. 

The standard type II migration: When the mass of 
planets in the range between Mgap < Mp < Merit (equa- 
tions ([28]) and ([29])). they undergo type II migration that 
proceeds as the gas disks evolve; the type II migration 
timescale Tmigji equals Tyis- Therefore, the planets move 
inwards with the velocity written as 

Vmigji — — — ■ (31) 

Slower type 11 migration: For planets with Mp > Merit i 
on the contrary, the type II migra tion rate slows down 
due to the inertia of the planets (|Sver &: Clarkd [19951 : 
llvanov et al.lll99"9[ ). That is why we refer this mode as 
to slower type II migration. As a result, the velocity of 
the planets becomes 

u 

Vmig,slowII — 7Z ' 7 «^ / 1. ^ 7: \'^^) 

+ Jmig,slowerIlMp/ jyicrit) 

where we have followed [Hellarv fc NelsonI (|2012D for tak- 
ing into account the effects of the inertia of planets. 
In addition, we have introduced a new free parameter 
fmig,siowerii- As showu below, both of the trapped type 

I and slower type II migration are important agents that 
regulate the radial distribution of planets in our model. 
Furthermore, it is currently uncertain how effective the 
inertia of planets is in slowing down the standard type 

II migration. Thus, it is useful to clarify the role of 
the slower type II migration by performing a parame- 
ter study wherein the value of fmig, slower ii varies (see § 
18. 3|) . We set f mig, slower ii = 1 for our fiducial model. 

6.3. Disk models 

We adopt the disk models discussed in §|31 More specif- 
ically, we use the characteristic surface densities at three 
disk inhomogeneities. In addition, the surface density of 
dust is required to examine planetary growth there. 



Table 4 

Values of fatg 

Dead zone Ice line Heat transition 

TFtf, 001 0.01 0.05 

We simply assume that 

S(J = fdtg'^g, (33) 

where f^ng is the dust-to-gas ratio. Table U) summarizes 
the values of fdtg at each disk inhomogeneities. We took 
fdtg = 0.01 at the dead zone, because the dust mass is 
canonically about a hu ndredth of the gas ma ss in pro- 
toplanetary disks (e.g. iDullemond et al.ll2Q07l ). At the 
ice line, condensation of water increases the dust density 
there and beyond. Therefore, we used fdtg = 0.05 at the 
heat transition. The reason that fdtg = 0.01 at the ice 
line is that we have already taken into account the effect 
of the ice line on by reducing the mean value of a (see 
equation (|T2|)). As a result, fdtg = 0.01 is reasonable for 
specifying there. 

6.4. Initial conditions 

We choose a value for the initial mass of cores 
~ O.OIMq, which is sufficiently smaller than the 
mass that is finally o btained by the oligarchic growth 
([Kokubo fc Idal [1991 [2002). We confirmed that this 
choice does not affect our results. 

The cores start growing at a position r at a time t. In 
principle, core formation takes place anywhere in disks. 
Nonetheless, we assume that the cores will quickly end 
up on one of the traps in the initial setup. It is noted 
that the assumption does not always assure the cores to 
be initially captured at their traps. This is because trap- 
ping happens only if the mass of the cores is larger than 
Mjnigj (see equation (|27pV Although one may consider 
the assumption of the initial r and r to be s omewhat 
artificia l, this is not the case. As shown by ilda fc LinI 
(|2008aD . planetary cores that undergo rapid type I mi- 
gration do not contribute to the population of gas giants, 
(since they plunge into their central star within the disk 
lifetime). This implies that only the cores that experi- 
ence slower type I migration will play an important role 
for reproducing the observed gas giants. Thus, it is rea- 
sonable to focus on planet formation proceeding only in 
planet traps that can substantially slow down the type I 
migration. 

Based on the assumption, it is only necessary to choose 
a distribution of the initial time r (or position r) for the 
growth of cores to begin. The positions of disk inho- 
mogeneities are related to the time r through the accre- 
tion rate (see equation (jTH])). Table [S] summarizes our 7 
choices of r which are selected to cover the entire disk life- 
time within which planetary growth and migration take 
place. For reference purpose, the initial positions that are 
determined by equations (0 , (0 , and ([8]) are also shown 
in the same table. It is noted that the multiple choices of 
the initial time result in forming multiple planets in each 
planet trap. We emphasize that the productivity of each 
planet trap - how many planets eventually form in each 
planet trap during the disk lifetime - and the relation 
between the number of finally formed planets and planet 
traps should be investigated separately. 
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We neglect the planet-planet interactions of the cores 
that grow in different planet traps - we leave this for our 
future work. 

6.5. Concurrent evolution of planetary growth and 
migration 

We may now follow the evolutionary tracks in the 
mass- semi- major axis diagram for planets that grow in 
all three planet traps. We adopt the above analytical 
prescriptions for planetary growth and migration. We 
summarize our technical procedures here. The stan- 
dard treatment of mass accretion and planetary growth 
is given in Appendix IB] 

When the mass of protoplanets is less than M^igj, 
their mass increases with time follo wing the standard 
oligarchic growth (see equation (jB2|) while their semi- 
major axes remain roughly the same. Time evolution also 
reduces the surface density of disks (gas and dust), and 
henc e th e growth rate also changes with time (see equa- 
tion (jB2p ). Once they acquire masses that are larger than 
M„iigj, they start to migrate inward. When they are at 
their planet traps, they move inward at the same rate as 
their traps. If they are left behind, they quickly catch up 
with their traps due to the standard rapid type I migra- 
tion, and then follow the movement of the traps. If their 
planet traps disappear due to convergence with other 
traps, it is assumed that the planets follow new planet 
traps that survive the convergence. If the planets become 
more massive than the critical mass of cores above which 
their envelope s ca nnot maintain hydrostatic equilibrium 
(see equation (|B7I) ). then accretion of gas onto the cores 
begins. The gas accretion ra tes a re regulated purely by 
the mass of cores (equation (jB8|) ). Through our exper- 
iments, we find that, for most cases, core formation is 
completed when they are captured in their traps. 

When a planet's mass reaches the gap opening mass, 
it undergoes the standard type II migration. This re- 
sults in " dropping-out" of the trapped planet from its 
planet trap and happens because the planet is now too 
massive to open up a gap in the disk (which leads to 
different orbital evolutions between the planet and the 
planet trap). When the planets are within the dead 
zone, type II migration becomes slower through a low 
value of a{— ao) while, for the planets outside the dead 
zone, the value of ua is used for the migration. If the 
planet attains the mass of fmax^max, where fmax is a 
controllable parameter (see Appendix iBt . its accretion is 
terminated. Mmax is a decreasing function of time (see 
equation so that planets that need a long time 

to grow up to gas giants tend to be less massive while 
planets that can quickly become gas giants tend to be 
more massive. Even when planet formation is largely 
complete, their disks may still have a sufficient amount 
of gas to drive type II migration. In this case, the type 
II migration is slowed down by the inertia of the planets. 
The accretion rate M declines with time, and at certain 
time M becomes equal to photoevaporation rates Mpe. 
When this is satisfied, the positions of the planets freeze 
in the mass-semi-major axis diagram. 

7. RESULTS 

We are now in the position to discuss the results of 
evolutionary tracks of planets that grow in disk inhomo- 
geneities. As shown in Fig. |4l most evolutionary tracks 



behave similarly. Therefore, we first discuss the results of 
the dead zone in detail (see § 17. and then examine all 
the three disk inhomogeneities (see § 17.21) . We compare 
the results with the observations in § 17.31 

7.1. Planetary growth in a dead zone trap 

The evolutionary track of a growing planet consists 
of four distinct phases in the mass-semi-major axis dia- 
gram (see Fig. [3]) . The first phase is formation of cores of 
gas giants through runaway and oligarchic grow th (e.g. 
IWetherill fc StewartI llQSfll : IKokubo fc 1<B [l998h . The 
mass of the core in this phase is high enough to keep 
up with the movement of the trap while the torque of 
the core acting on the disk is too weak to open up a gap 
there. Thus, the core remains within a trapping regime 
and follows its movement. The timescale of this phase is 
order of ^ 10^ — 10^ years, which is much shorter than 
the disk lifetime (r^isfc ~ 8. 8 x 10^ years in this setup) as 
shown in previous studies (|Kokubo fc Idall2002D . Hence, 
the protoplanet moves upwards in mass while moving 
little in orbital radius or period. 

As the feeding zone empties, core formation is termi- 
nated and the second phase begins, wherein gas accretion 
onto the envelope occurs. It was well known that the 
timescale of this phase w as problematically l ong for ear- 
lier models (> 10'' vears. iPoUack et all 119960 . However, 
recent studies improved the previous models and revealed 
that the timescale is highly sensitive to the optical depth 
of the envelope. For these realistic condi tions, it is sig- 
nificantly shorter than the disk lifetime (jLissauer et al.l 
120091) . In our calculation, this timescale is about 2 x 10 
years and hence the core still has sufficient time to finally 
grow up to a gas giant within Tdisk ■ The core during most 
of this phase is trapped. As a result, its radial evolution 
is mainly determined by the slow movement of the dead 
zone trap, and the protoplanet moves to shorter radii and 
periods while at nearly a constant mass. Toward the end 
of this phase, the core becomes massive enough, so that 
the tidal torque it exerts upon the disk becomes compa- 
rable to the viscous torque that evolves gas disks, leading 
to gap formation in the disks and type II migration of 
the planetary core. 

When the mass of the gaseous envelope cannot be sup- 
ported by the gas pressure, runaway gas accretion onto 
the core takes place (Phase HI). The timescale of this 
phase is very short (< lO'^ years), and consequently its 
evolutionary path is almost vertical in the mass-semi- 
major axis diagram. These three successive phases are 
the main path to forming gas giants in the core acc retion 
scenario ((Pollack et al.lll996t iLissauer et al.ll2009D . The 
massive planet opens up a gap in the disk and undergoes 
type II migration. This switch from type I to type II 
migration results in "dropping-out" of the planet from 
the moving trap and decouples it from the movement of 
the planet trap. 

The onset of Phase IV completes the formation of a gas 
giant. During this phase (> 10^ years), type II migra- 
tion moves the gas giant inward further. However, this 
process is minimized by the inertia of th e massive planet 
(|Sver fc ClaAa[T99l llvanov et al.l[T999l) . 

Planets arrive at their final position in the mass- 
semi-major axis diagram when the disk is finally dis- 
sipated. Photoevaporation of the disk by high en- 
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Table 5 

The initial times and positions 



Tiic initial time (yr) 


Dead zone (AU) 


Ice line (AU) 


Heat transition (AU) 


10^ 


7.3 


11.7 


42.4 


2 X 10^ 


5.7 


7.4 


22.3 


4 X 10^ 


4.4 


4.6 


11.7 


8 X 10^ 


3.4 


N/A 


6.1 


1.6 X 10*^ 


2.6 


N/A 


3.2 


3.2 X 10*^ 


2.0 


1.9 


N/A 


6.4 X 10^ 


1.6 


1.2 


N/A 



10^ 



10" 



10"' 



10-' 



I I I 1 1 1 1 1 1 — I I I 1 1 1 1 1 1 — I I I 1 1 1 1 1 1 — I — I I 1 1 1 



Phase IV 

9^ 



Phase 



We assume that planet formation does not take place in a planet trap when the planet trap disappears due to 
convergence with a dead zone trap. N/A represents such cases. 

ergy radiation from the central star is likely to be the 
dominant mechanism of gas dispersal in the disks (e.g. 
IGorti fc^oUenbach 2009.) , and will terminate type II mi- 
gration. As a result, the final orbital period and mass of 
the planet are achieved. Thus, Fig. [3] summarizes how 
concurrent evolution of planetary growth and migration 
proceeds in the mass-semi-major axis diagram: a core is 
formed in a dead zone trap that is initially located at 
~ 7 AU. Following the movement of the dead zone trap, 
the core is transported to ~ 3 AU. Simultaneously, it 
undergoes the two main phases of gas giant formation. 
The completion of the final runaway gas accretion onto 
the core and subsequent type II migration involve further 
evolution of the planet in the diagram. When photoevap- 
oration becomes important, the gas disk is removed and 
the position of the planet in the diagram " freezes-out" . 



Phase 




Phase I 



- Dead zone 

I I I I 



_i I I I I I I I I I 



0.01 



0.1 



10 



100 



7.2. Planetary growth in all the three planet traps 

Fig. |4]shows the computed evolutionary tracks of plan- 
ets that grow at all three disk inhomogeneities. Different 
lines at each planet trap correspond to different evolu- 
tionary tracks in which planetary growth starts at dif- 
ferent times (see Table [S]) . Despite the difference in the 
starting time (and position) , most planets formed at the 
dead zone and heat transition traps end up at r 1 AU 
500 days) and r ^ 0.1 AU (~ 10 days), respectively. 
At the heat transition trap, the surface density of dust is 
low. Therefore, cores that grow there spend a long time 
in the trapping phases (Phase I and II). This maximizes 
the distance over which cores are transported and re- 
sults in the distribution of cores that hover preferentially 
around > 1 AU. Since the low mass cores get distributed 
over smaller orbital radii and less time remains for the 
cores to grow up to gas giants, they finally remain less 
massive (< lOOM®), and are located around smaller or- 
bital radii 0.1 AU). The same argument is applied 
to planets formed in the dead zone trap. However, the 
surface density of dust at the dead zone is considerably 
higher than that at the heat transition. Consequently, 
the final mass of cores trapped at dead zones becomes 
larger, core formation completes earlier, and the distri- 
bution of cores is shifted to ~ 3 AU. These combined dif- 
ferences result in the populations of more massive planets 
orbiting at ^ 1 AU. 

The evolutionary tracks associated with protoplanets 
carried by the ice line trap show some differences. The 
resultant planetary population spreads out over a wider 
range in the mass-semi- major axis diagram (see Fig. 2]). 
Nonetheless, this can be also understood by the the sur- 
face density of dust and the resultant core formation 



r (AU) 



Figure 3. An evolutionary track of a planet that grows in a dead 
zone trap. The track (denoted by the thick line) can be divided into 
four phases. In Phase I, core formation takes place very rapidly in 
T ~ 10^ — lO" years, which is much faster than the radial movement 
of the trap at that time. This results in largely vertical motion in 
the diagram. In Phase II, the core accretes gas onto its envelope. 
Its timescale is very slow (~ 2 X 10® years in this case). Therefore, 
it moves horizontally in this diagram. The mass of the core in 
Phase I and most of Phase II is within the trapping regime that 
is represented by the upper and lower thin lines. Toward the end 
of Phase II, the core drops-out from the trap by opening up a 
gap in the disk and undergoing type II migration. Phase III is 
runaway gas accretion onto the core. The timescale of this phase 
is very short (< 10^ years). As a result, it moves vertically in 
this diagram. Planet formation completes during Phase I to III. 
In Phase IV (> 10® years), the gas giant moves inward due to 
type II migration that is slowed down by the inertia of the planet. 
When photoevaporation of the gas disk becomes important, type 
II migration is terminated and its final radial position and orbital 
period are obtained. 



there. At the ice line, the surface densities are sub- 
stantially higher than that at the dead zone and heat 
transition and hence the formation of cores is most effi- 
cient. This typically results in most massive cores. At the 
early stage of disk evolution, therefore, the most massive 
cores are preferentially formed there. They can readily 
drop-out from the moving trap and pile up around larger 
orbital radii (r ~ 5 AU). These massive cores at larger 
orbital radii lead to the formation of more massive gas 
giants that finally orbit at > 1 AU. In the later stage of 
disk evolution, the high dust densities at the ice line can 
still form cores while at that time the other traps not due 
to lower dust density there. This is the physical reason 
of the wide spread of planetary population due to the ice 
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Figure 4. Evolutionary tracks of planets that grow in all three 
planet traps. The tracks for the dead zone are denoted by the 
red thick lines, the ice line by the green, and the heat transition 
by the light-blue. Corresponding thin lines represent the trapping 
regimes. Different tracks correspond to planetary growth that ini- 
tiates at different times (see Table O . The transport mechanism 
of cores by planet traps plays the crucial role in producing the 
mass-period relation; low mass cores that need longer time to grow 
are more likely to be transported toward smaller orbital radii while 
massive cores that can readily drop out of the moving traps tend 
to distribute further away from the star. Thus, there are distinct 
populations that arise from the difference in the properties of the 
planet traps and the resultant planetary growth, which results in 
the t rend that planetary ma ss increases with period. Earlier pa- 
pers, |ld^^XI3 l|2004l . [2008bl ) predicted a planet desert demarcated 
by the black rectangle. We emphasize that our model predicts the 
presence of planets there. 

line traps. 

7.3. Comparisons with the observations 

We now compare our results with the observations. As 
already presented in Fig. |4l our model shows that the 
superposition of all tracks for planets that grow in three 
planet traps constitutes a theoretical mass-period rela- 
tion, wherein the final distribution of the mass of the 
planets is an increasing function of their periods. This is 
consistent with the observed mass-period relation, as the 
observational data scatter around the locus of end points 
of our tracks (see Fig. [5|). 

This is one of the most important findings in this pa- 
per. As discussed in § 17.21 this arises from the fact that 
there are considerable differences in the properties of the 
planet traps that regulate planet formation and migra- 
tion. As a result, different planet traps have different 
preferred loci at which evolutionary tracks end up in the 
mass-semi-major axis diagram. Thus, planet traps act as 
a filter for distributing cores - massive cores readily drop 
out from moving traps and tend to orbit further away 
from the central star while low-mass cores are trapped 
for a long time and tend to orbit close to the host star 
- and play the central role in generating the theoretical 
mass-period relation. 

In addition, the prediction that distinct sub- 
populations can arise depending on the trapping mecha- 
nism has several observational consequences. For exam- 
ple, our model provides a physical explanation for the 
observed pile up of gas giants at ~ 1 AU. This again 



relies on the argument that planet formation efficiency 
highly depends on the surface density of dust at planet 
traps. At the dead zone and ice lines, the dust density 
is expected to be high due to the low disk turbulence, 
and hence planet formation rates are high there. On the 
other hand, the formation rate would be low at the heat 
transition trap due to low dust density. This results in 
a general trend that more planets are readily formed at 
the dead zone and ice line traps that end up at r ~ 1 AU 
(see Fig. [5]). 

Furthermore, our model predicts the population of low 
mass planets (< SOAf®) with r < 0.5 AU. This arises 
from planet formation that takes place in the moving 
ice line trap (see Fig. [5]). Even in the later stage of 
disk evolution, the highest dust density there enables the 
formation of low-mass planets that end up in the desert. 
On the contrary, the most advanced pop ulation synthesis 
models predict a planet desert there (jida fc LinI 12004 
I2008bl also see the footnote 3 in §[!]). The presence of 
the many observed exoplanets in the region agrees well 
with our findings. 

Finally, our models predict the existence of planet 
deserts that are quite different in the mass-peri od space 
than those claimed bv llda fc Lii] (pOOl [2008H ) . Fig. M 
shows our deserts, denoted by hatched regions. They 
are produced due to trapping and subsequent transport 
of cores. This leads to the evacuation of the cores from 
these regions in which they have initially grown up. As 
a result, these regions are regarded as void of planets. 
More specifically, we define our deserts by estimating 
the mass ranges of planets that can be captured at the 
planet traps and following their movement: Mp < Mgap 
and Tmig.i < Tyis (see §[5]). This kind of planet desert is 
active only for gas disks. There are a number of possi- 
bilities to fill out our deserts: that successive formation 
of rocky planets after gas disks disperse may ultimately 
fill out the regime; that, even in the epoch of gas disks, 
planetary cores formed far beyond our deserts may even- 
tually distribute there due to planetary migration; and 
that planet-planet scatterings induced by convergence of 
multiple planet traps may deliver the scattered cores into 
our deserts. Nonetheless, our predictions are valuable in 
a sense that such regions are the primary target of the 
current and ongoing obse rvational surveys (jMavor et al.l 
[Ml [Howard et al.li2012l ). 

8. PARAMETER STUDIES 

We perform parameter studies by varying disk and stel- 
lar parameters in order to examine how robust our find- 
ings discussed in § 17.31 are. Also, we investigate the ef- 
fects of the inertia of planets by changing the value of 
fmig,siowerii in § 18.31 to differentiate them from the role 
of planet traps discussed above. 

8.1. Disk parameters 

We first focus on parameters for dead zones. We have 
adopted the parameterized treatment for the structures 
of dead zones, wherein they are represented by Sao and 
SA (see equation ([3])). Even in the most recent stud- 
ies, it is still somewhat unc ertain what the precise struc- 
ture of the dead z ones is (jMatsumura fc Pudritd [20061 : 
[Martin et an[2012D . Therefore, we utilize our parameter 
study in order to discuss how sensitive our findings are 
to the structures of the dead zones. 
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Figure 5. Compariso ns with the observat ions. The observational 
data are adopted from lCTavor et all I I2011I) (as Fig. [TJ. Our the- 
oretical mass-period relation is consistent with the observations. 
Also, the presence of many observed low mass planets (< SOM^) 
at r < 0.5 AU provides further support on our model. 
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Figure 6. Prediction of planet deserts for the CTTS case. The 
desert produced by the dead zone trap is denoted by the red 
hatched region, the ice line trap by the green, and the heat transi- 
tion by the light-blue. These regions are the consequence of trap- 
ping of type I migr ators and defined by the gap opening mass Mgap 
(see equation I I28I I) and the mass of planets above which they can 
keep up with the movement of their planet traps (see equation I I27I I , 
also see Table [Sjl. 



Table |6] summarizes parameters we varied. For Runs 
Al and A2, the value of T,ao is changed while sa varies 
for Runs A3 and A4. Any other parameters remain the 
same as the fiducial ones for all the four runs. Fig. [7] 
shows the results of the evolution of the positions of three 
disk inhomogeneities (the left column), the evolutionary 
tracks of planets (the central column), and the trapping 
regimes (the right column). The top panels are for the 
case of Run Al, the second for the Run A2, the third 
for Run 3, and the bottom for Run4. One immediately 
observes that the results for all the four cases, especially 
the behaviors of the evolutionary tracks, are surprisingly 
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similar to those of the fiducial case and give all three 
key results. Thus, we can conclude that our findings 
discussed in § 17.31 are robust even if the structures of 
dead zones somewhat change due to their surrounding 
environments and disk configurations. 

8.2. Effects of Stellar mass 

The above parameter study leads to the conclusion 
that disk inhomogeneities and the resultant multiple 
planet traps play the crucial role in reproducing the key 
properties of observed exoplanets. Nonetheless, there are 
few populations of planets that are not covered by the 
fiducial model: gas giants orbiting at r > 5 AU. We now 
examine whether or not this population is also predicted 
by our model. In order to proceed, we change the stel- 
lar mass from 0.5 to 0.9Mq and keep other parameters 
the same as the fiducial ones. The main motivation for 
changing the stellar mass is that the observational data 
arc obtained from low- to high-mass stars that cover F, 
G, and K stars. 

Fig. [8] shows the results of the movement of the disk 
inhomogeneities, the evolutionary tracks of planets that 
are formed in the traps, and the locations of the planet 
deserts. This figure confirms that the planets not covered 
in our previous setup (lO^M^ < Mp < 5 x lO^M^ and 5 
'^r < 10 AU) can indeed be reproduced. This is be- 
cause high-mass stars result in high accretion rates (see 
equation (|19|) ). which corresponds to the situation that 
disks have high mass. As a result, planet formation effi- 
ciencies at all three planet traps become high, and most 
formed planets readily attain high mass, which ends up 
with planets distributing further away from the central 
star. Thus, this finding indicates that the full range of 
the statistical properties of exoplanets can be understood 
by our model. 

8.3. Role of slower type II migration 

Based on the above discussion, a number of the fun- 
damental statistical properties of observed exoplanets 
are most likely to be explained by multiple planet traps 
that capture and transport planetary cores, depending 
on planetary mass. As shown in Fig. |4l however, slower 
type II migration also drives radial drifts in the mass- 
semi-major axis diagram, following evolutionary tracks of 
growing planets^ In addition, the slower type II migra- 
tion is a function of planetary mass (see equation (p2|) V 
Thus, they imply that a combination of planet traps and 
slower type II migration (not only planet traps) may play 

* The standard type II migration also plays some role in changing 
the radial distribution of planets. However, its effect is minimal, 
because runaway gas accretion proceeds so rapidly that, once plan- 
ets obtain the gap-opening mass, they immediately achieve their 
mass above which the inertia of the planets is effective. As a result, 
the time interval for the standard type II migration is short enough 
to neglect the effects. 
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Figure 7. Parameter study of dead zones (see Table (6)1. The movements of three planet traps are shown in the left column (as 
Fig. [2J, the evolutionary tracks of planets that grow there in the central (as Fig. O, and the behaviors of the trapping regimes in 
the right (as Fig. [BJ. The top panels show the results of Run Al, the second for Run A2, the third for Run A3, and the bottom 
for Run A4. The results are quite similar to those of the fiducial model, and hence our important findings such as origins of the 
observed mass-period relation and the pile up at 1 AU and a prediction of low-mass planets with tight orbits are maintained for 
dead zones that can have a variety of structures. 



a central role in reproducing the observations. In order 
to examine the effects of the slower type II migration on 
our results, we carry out a parameter study in which the 
value of fmig.siowerii chaugcs (see Table [71). 

Fig. [5] shows the results of the evolutionary tracks of 
planets growing in all the three planet traps. The top 
left panel shows the results of Run Bl, the top right for 
Run B2, the bottom left for Run B3, and the bottom 
right for Run B4. For two top panels, there is no mass 
dependency in equation ((32l) while it depends on plane- 
tary mass for two bottom panels. Comparison of Run Bl 
(Top left) with any other runs demonstrates that some 
kind of process which slows down the standard type II 



migration is clearly needed for reproducing the observed 
population of exoplanets even if planetary cores are saved 
due to planet traps. This is because the standard type II 
migration that proceeds as local viscous timescale leads 
planets to spiraling into the host stars within the disk 
hfetime at r < 10 AU. 

When the type II migration somewhat slows down 
sufficiently, our findings discussed above, especially the 
mass-period relation, are valid for various cases (see 
Top right and two Bottom panels). It is important 
that the overall feature of the resultant populations is 
very insensitive to the origins of slowing down mecha- 
nisms (mass dependence vs independence, compare top 
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Figure 8. Parameter study of the variation of stellar mass. The movements of three planet traps are shown in the left panel, the 
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right panel with Fig. [5]), ahhough the radial distri- 
bution somewhat varies for each case. We confirmed, 
through experiments, that the observed planetary pop- 
ulation can be reproduced when fmig.siowerii ^ 1 or 

fmig,slowerII — fMcrit/Mp with / > 1. 

In addition to the slower type II migration, pho- 
toevaporation of gas is also important for estab- 
lishing the final radial distribution of planets, since 
it terminates the migration. The combined ef- 
fects of type II migration and photoevaporation 
on planetary populati ons were dis c ussed in the lit- 
erature (e.g. Armitag e et al.l [20021: iMatsuvama et alj 
200a Armitagc 200 71: [Alexander fc Armitaed l200a 
Alexander &: Pascucciil2012f) . Our models are more fun- 



damental than theirs in a sense that we have incorporated 
the effects of type I migration which can be captured at 
planet traps. We will investigate in detail the combined 
effects in the forthcoming paper. 

In summary, we can conclude that planet traps, not 
slower type II migration, play the primary role in repro- 
ducing the observations, provided that the standard type 
II migration slows down sufficiently when planetary mass 
exceeds the local disk mass. 

9. DISCUSSION & CONCLUSIONS 

We have constructed and followed the evolutionary 
tracks of protoplanets generated by combining core ac- 
cretion together with the movement of planet traps in 
evolving disks. We have focused on three types of in- 
homogeneities in protoplanetary disks and the resultant 
planet traps: dead zones, ice lines, and heat transitions. 

We have demonstrated that the plant traps play two 
fundamental roles in planet formation and migration. 
The first is to trap cores of gas giants that otherwise 
undergo rapid type I migration. This trapping leads to 
the formation of gas giants orbiting at 0.01 AU ^ ^ 10 
AU without the cores falling into the central star before 
disk evolution is terminated by disk photoevaporation. 
The second is to transport the trapped cores very slowly 



from large to small periods. This transport distance is 
regulated by the gap opening mass and the properties 
of the planet traps and hence is well coupled with plan- 
etary growth histories; different planet traps result in 
different efficiencies of planet formation and migration. 
Consequently, planet traps are regarded effectively as a 
filter for distributing cores - massive cores tend to hover 
around large periods while low-mass cores around short 
periods. 

We have seen that the combination of planet traps, 
planetary growth and type II migration in evolving disks 
generates planetary populations that have a wide range 
in mass and period. The final positions of planets are 
determined when photoevaporation of gas disks (rather 
than disk viscosity) plays the dominant role in disk evo- 
lution. It is noted that the final distribution of planets 
can be largely affected by the combination of slower type 
II migration and the photoevaporation rate that is ad- 
justed by fpe (see equation pO| and Tablc[2]). In order to 
examine this dependency, we will perform a more com- 
prehensive parameter study in a subsequent paper. In 
addition, different accretion histories of planets will re- 
sult in differences in the composition of the planets and 
their atmospheres, which can be investigated by extend- 
ing our models. 

What about purely dynamical effects arising from 
planet-planet interactions? These are well known to 
be important for unde rstanding the observ ed eccentric- 
ity distribution ( e.g. iRasio fc FordI I1996D . However, 
iMatsumura et al.l (|2010l )" have recently clarified through 
numerical simulations that the semi-major axes of plan- 
ets are determined mainly by planetary migration in the 
gas disks while the eccentricities are determined after the 
gas disks dissipate severely. This indicates that the re- 
sults of planetary evolution in the gas disk phase will not 
be washed out by subsequent planetary dynamics. 

We list our major findings below. 

1. We have demonstrated that the wide range of end 
points of planets evolving in the mass-semi-major 
axis diagram establishes a theoretical mass-period 
relation in which planetary mass is an increasing 
function of orbital period (see Fig. [S]). This is in 
excellent agreement with the observational data in 
a sense that the data scatter around the end point 
of our evolutionary tracks. 

2. We have shown that the many tracks of dead zones 
and ice lines preferentially tend to end up at r^ 1 
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Figure 9. Parameter study of the effects of slower type II migration (see Table O. The resultant evolutionary tracks of planets 
that grow in all the three planet traps are shown (as Fig. (SJ. The top left panel denotes the results of Run Bl, the top right for 
Run 132, the bottom left for Run B3, and the bottom right for Run B4. For two top panels, the slowing rate is independent of 
planetary mass while it depends for two bottom panels. When no slowing mechanism is in action (Run Bl, fmig,slowerll = 0)> 
the observed population of exoplanets cannot be reproduced even if planetary cores are captured and distribute by moving planet 
traps. The observations can be reproduced only if the standard type II migration sufficiently slows down {fmig,3lowerll 1 or 
fmig, slower II = f^'Icrit/^p with / > 1). It is important that mass dependency on the slowing mechanism of the type II migration 
does not matter. Thus, this suggests that planet traps are more crucial for understanding the statistics of observed exoplanets. 



AU (see Fig. [5]). Combined with an argument that 
planet formation efhciencies are reasonably high 
there, the preference provide a physical explana- 
tion for the pile up of observed gas giants at ^ 1 
AU. 

3. We have also demonstrated that planets that grow 
in dead zone traps end up at r ~ 1 AU, ice line 
traps at 0.03AU ^ r < 3 AU, and heat transition 
traps at r ~ 0.1 AU (see Fig. S]). The resulting 
wide range of planets in the mass-semi-major axis 
diagram is insensitive to the detailed structure of 
dead zones and accounts for a number of important 
observational trends. 

4. We have also shown that moving ice line traps can 
put planets in the planet deserts that were pre- 
dicted by the earlier population synthesis models 
(see Fig. [5]). As denoted by Fig. [U the desert is 
located in the range of planetary masses (5-50 M0) 
and semi-major axes (0.04-0.5 AU). The recent ob- 
servations discover many planets in the deserts. 
Thus, our models are more consistent with the ob- 
servations in this regard. 



5. We predict planet deserts that arise from the na- 
ture of planet traps. They have physically different 
origins from the earlier claimed ones. Our deserts 
are relevant only when protoplanetary disks have 
sufficient amount of gas that drives rapid type I mi- 
gration. This suggests that the deserts can be filled 
by planets that form far beyond the deserts and 
eventually migrate there in gas disks. The planets 
are likely to emerge after the gas disks severely dis- 
perse and/or to be transported there due to planet- 
planet scatterings. Our deserts are present in the 
range of planetary masses (1-50 M^) and semi- 
major axes (1-10 AU), which covers the primary 
target of ongoing and future observational surveys 
such as the HARPS and Kepler missions. 

6. The more massive the host star, the more the evo- 
lutionary tracks in the mass-period diagram are 
pushed towards large disk radii. This arises be- 
cause of the much more rapid accretion rates in 
more massive systems (Af cx Af^). 

In the forthcoming paper, we will use N-body simula- 
tions to take into account the physics of planet-planet 
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interactions that can be induced by growing planets in 
different planet traps. 
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APPENDIX 

A: CHARACTERISTIC MASSES IN THE MASS-SEMI-MAJOR AXIS DIAGRAM 

We discuss the segmentation of a diagram for planetary mass versa semi-major axis. As an example, the bottom 
panel of Fig. [10] (Left) shows evolution of four characteristic masses (Mmax, Merit, Mgap, and Mmigj) at a dead zone 
for the fiducial case (also see Table [3|). Every position of the dead zone that is specified by the time t defines four 
masses (see circles on Fig. [TU] (Left) as an example). As the dead zone moves inwards following the disk evolution 
(see the top panel of Fig. [TU] (Left)), these four masses also move inwards at the same rate. As a result, four lines are 
drawn that track the evolution of the characteristic masses in the diagram. The top line denotes Mmax, the second 
for Merit, the third for Mgap, and the bottom for Mmigj- Thus, the top and third lines define the boundaries of type 
II migration in the mass-semi major axis diagram while the third and bottom one defines the boundaries for trapping 
type I migration. The second line defines where the inertia of planetary mass becomes effective that can slow down 
type II migration. 

The bottom panel of Fig. [TU] (Right) shows the evolution of four characteristic masses at all three disk inhomogeneities 
for the fiducial case. The type II regimes for each inhomogeneity are denoted by the coarse hatch regions while the 
fine hatch ones represent the type I trap regimes. The solid lines denote Merit- The dead zone is denoted by red, the 
ice line by green, and the heat transition by light blue. Disappearance and re-appearance of these regions correspond 
to the behavior of each inhomogeneity that is shown in the top panel. 

B: ANALYTICAL PRESCRIPTIONS FOR PLANETARY GROWTH 

We summarize the standard results of the core accretion scenario needed to follow the growth of planets as they 
move in disks. 



Stage I: Formation of cores 

Formation of cores in planetesimal disks is well understood in the literature (e.g. iKokubo &: Idal[2002| ). and the 
growth timescale in the disks is given as 
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where S^; is the surface density of dust. Me is the mass of a core, 6 = 10 is a parameter for determining the feeding 
zone of the core (see below), and m = lO^^g is the mass of planetesimals that are accreted onto the core. Adopting 
this timescale, the growth of cores is regulated by 



dMp 
dt 



(B2) 



Stage II: slow gas accretion onto the envelopes 

As cores grow, their feeding zones Arc empty. These zones are scaled by brn, where rn = (Mc/(3M,))i/'^ is the 
Hill radius of a core and b ^ 10. The decrease of planetesimals in the zones results in the reduction of the accretion 
rate of cores Me- In the limit of the modest to high velocity dispersion a of planetesimals that are accreted onto the 
cores, Me can be written as (|Safronovlll972l : llda fc Lin. 20041 
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where Re is the radius of cores. Planetesimals within the feeding zones can reach cores when tr/fi ^ Are{— brn) 
Assuming Re to be similar to that of the Earth; 



Re 



6.4 X 10* cm 



Me 

M(. 



1/3 



5.5 g cm 



-3 



-1/3 



(B4) 
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Figure 10. Left: (Top) The movement of a dead zone trap in evolving disks. The initial position is determined by assuming a 
start time of 10^ years. The final position arises when the accretion rate M equals the photoevaporation rate Mpe, - a condition 
which also defines the disk lifetime {rdisk — 8.8 X 10® years in our fiducial model). (Bottom) The evolution of four characteristic 
masses at the dead zone trap. The top line represents Mmax, the second for Merit, the third for Mgap, and the lowest one for 
Mmigj (see TablelSj. Every position of the dead zone (shown in the top panel) defines four masses (see some symbols as examples), 
and as the dead zone moves inward, they also move inward. Right: (Top) The movement of all three traps in evolving disks. The 
dead zone trap is denoted by red, the ice line by green, and the heat transition by light blue. They move inward due to viscous 
evolution of the disks, but at different rates. (Bottom) The evolution of four characteristic masses at all three traps. The color 
scheme is the same as the top panel. The type II regimes that are defined by Mmax and Mgap, are denoted by the coarse hatched 
regions while the trapping regimes that are specified by Mgap and Mmigj are represented by the fine hatched regions. The solid 
lines denote Merit- The disappearance and re-emergence of these regimes correspond to those of the positions of the planet traps 
that are shown in the top panel. 

Mc is given as 



M, - 3.0 X 10-«Me yr"! ( ^) 3 

'10/ \5.5 g cm 



-1/3 ^MrV fAL^-'" 



Ma 



Ma 



10 g cm 



When all the planetesimals in their feed ing zones are consumed, the cores att ain the maximum mass that is known 
as the isolation mass, which is defined by (IKokubo fc Idal[200l llda fc Linl[20Ql 



Mr,.,so = 27rr A r.Srf ~ O.WMa, I — 



3/2 



^ \ 3/2 / \ 3 / ,^ \ -1/2 



10 g cm 



M, 



(B6) 



The accretion of gas onto cores and subsequent envelope formation are initiated when the mass of core becomes larger 
than 

M,,,„, 10/,,e„tMe Lo-6Mtyr-0 ■ (^^^ 

This critical mass was originally derived from a series of numerical si mulations that inve stigated the effect of cores' 
accretion rates and opacit y in the e nvelop e on formation of gas giants (jlkoma et al.ll2000l ). Here, we have adopted a 
simplified one, followin g Ida fc LinI (|2004D . Recent studies, however, re vealed that Mc^crit might be smaller than that 
predicted by equation (IB7p with fc,crit = 1 (e.g. iHori fc Ikomal[2"011[ ). In order to take this into account, we have 
introduced a dimensionless factor /c crit- 
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Table 8 








Important parameters for planetary growth 






Symbols 


Meaning 




Value 


fc,crit 


A factor linked to the critical mass of cores (equation 


l|B7ll) 


0.3 


c 


Exponent of the Kelvin-Helmholtz timescale (equation 




8 


d 


Exponent of the Kelvin-Helmholtz timescale (equation 


m) 


2.5 


fmax 


A factor linked to the maximum mass of planets 




0.1 



The gas accretion rate of cores is prescribed by (|Ida fc Liijl200l 



dMp Mp 



dt TKH 

where the Kelvin-HehTiholtz timescale is given as 



(B8) 



-d 



where c = 8 and d = 2.5. This is a simphfied timescale that was origina lly estimated from numerical simulation s 
(jlkoma et al.ll2000l ). As shown by the more detailed numerical simulations (jPoUack et al.l[l996l : iLissauer et aLll2009( ). 
this stage is slow (> 10^ years). 

Stage III: runaway gas accretion onto the cores 

When planets become massive enough, runaway gas accretion onto the cores starts. In the detai led numerical 
simulation s, this s tage commences when the envelope of cores becomes more massive than the cores (jPoUack et al.l 
[1996: .Ikoma et al.ir2000: .Lissauer et al...2009i) . Nonetheless, we monitor this stage by the condition that 

^"^^ < 1. (BIO) 



105 yr 



This is because our approach i s ra ther simple. This condition never affects our results. The growth rate of this stage 
is also prescribed by equation (jB8[) . 

It is totally unclear how gas accretion onto the cores terminates and what physical process(es) determines the final 
mass of gas giants. Therefore, we assume that Stage III continues until planets gain the mass fmaxMmax, where fmax 
is an adjustable parameter. 



Parameters for planetary growth 

As discussed above, planetary growth and consequent evolutionary tracks of planets are regulated by four parameters 
in our model (see Table [8]). The parameter fc,crit governs the onset of gas accretion of cores, a set of parameters c and 
d determine the efhciency of gas accretion onto cores, and fmax contro ls the final mass of planets. We denote these 
values given in Table [5] our fiducial model. Note that llda fc LinI ()2004[ ) adopted the values of fc,crit = 1, c = 9, and 
c? = 3 rather than our set of fc,crit — 0.3, c = 8, and d = 2.5. We performed a parameter study and confirmed that 
different choice of these values does not change our results v ery much. There f ore, ou r choice and results are robust 
in a sense that we can compare our results with the results of llda fc LinI (|2004L [2b08bf ). We present only a parameter 
study of fmax in Appendix [Cl because the choice of this value is probably the most uncertain. 

C: A PARAMETER STUDY FOR PLANET GROWTH 

As discussed in Appendix [B] our choice of most parameters that regulate planetary growth is based on the physical 
consider ations and the m ore recent results of detailed simulations. Hence our choice is compatible with the original 
work of (|Ida &: Linll2004f ). However, there is one exception, which is fmax which constrains the maximum mass of 
planets. This is mainly because there are no firm physical arguments and simulations of how the final mass of planets 
is established. Therefore, we now examine how the value of fmax affects our results. 

Table [9] summarizes our parameter study on fmax- For Run Cl, we took fmax = 0.03 while fmax = 0.3 for Run 
C2. Except for the value of fmax, we adopted the same values of the fiducial model. Fig. [TT] shows the results of the 
evolutionary tracks of planets for both cases. The left panel shows the results of Run Cl while the right one for the 
Run C2. The results of both cases are generally very similar to that of the fiducial model. More specifically, both cases 
produced theoretical mass-period relations that are broadly consistent with the observations. On closer examination, 
we see that the pile up of gas giants at ^ 1 AU and the presence of low-mass planets with small orbital radii are 
relatively affected by the value of fmax- This is indeed expected. If fmax has a small value, then low mass planets 
become the main product and planet formation completes earlier. Consequently, the planets experience substantial 
inward type II migration. In addition, the effect of the inertia of the planets that slows down the type II migration 
is also reduced. This results in larger populations of low-mass planets at small orbital radii. Thus, models with small 
values of fmax have difiiculty in reproducing the pile up at 1 AU. The opposite happens for large values of fmax- In 
this case, massive planets are preferentially formed and the completion of planet formation takes place at later time. 
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Table 9 

Parameter study of planetary growth 

fmax 

Run CI 0.03 
Run C2 0.3 
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Figure 11. Parameter study of the variation of maximum mass of planets fmax (see Table|9]|. The evolutionary tracks of planets 
that grow in all three planet traps are shown (as Fig. [Sj. The result of Run CI (fmax = 0.03; low mass case) is presented in the 
left panel while that of Run C2 {fmax = 0.3; high mass case) in the right. Although the pile up at 1 AU and low mass planets with 
small orbital radii are slightly affected by fmax, the resultant populations are still well understood by the physical considerations 
discussed in §[3 

SO that the subsequent inward type II migration is significantly suppressed due to both the shorter remaining time 
and the larger inertia of planets. As a result, the population of low-mass planets with tight orbits declines while the 
1 AU pile up is enhanced. 

In summary, the results depend slightly on some basic parameters such as fmax- Nonetheless, they are well under- 
stood by the physical arguments presented in § [71 Therefore, our findings are reasonably robust for a wide range of 
the parameter space. 
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